Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture.

BACKGROUND: The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function. RESULTS: We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection. CONCLUSIONS: We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Complex Feline Disease Mapping Using a Dense Genotyping Array.

The current feline genotyping array of 63 k single nucleotide polymorphisms has proven its utility for mapping within breeds, and its use has led to the identification of variants associated with Mendelian traits in purebred cats. However, compared to single gene disorders, association studies of complex diseases, especially with the inclusion of random bred cats with relatively low linkage disequilibrium, require a denser genotyping array and an increased sample size to provide statistically significant associations. Here, we undertook a multi-breed study of 1,122 cats, most of which were admitted and phenotyped for nine common complex feline diseases at the Cornell University Hospital for Animals. Using a proprietary 340 k single nucleotide polymorphism mapping array, we identified significant genome-wide associations with hyperthyroidism, diabetes mellitus, and eosinophilic keratoconjunctivitis. These results provide genomic locations for variant discovery and candidate gene screening for these important complex feline diseases, which are relevant not only to feline health, but also to the development of disease models for comparative studies.

Early onset adult deafness in the Rhodesian Ridgeback dog is associated with an in-frame deletion in the EPS8L2 gene.

Domestic dogs exhibit diverse types of both congenital and non-congenital hearing losses. Rhodesian Ridgebacks can suffer from a progressive hearing loss in the early stage of their life, a condition known as early onset adult deafness (EOAD), where they lose their hearing ability within 1-2 years after birth. In order to investigate the genetic basis of this hereditary hearing disorder, we performed a genome-wide association study (GWAS) by using a sample of 23 affected and 162 control Rhodesian Ridgebacks. We identified a genomic region on canine chromosome 18 (CFA18) that is strongly associated with EOAD, and our subsequent targeted Sanger sequencing analysis identified a 12-bp inframe deletion in EPS8L2 (CFA18:25,868,739-25,868,751 in the UMICH_Zoey_3.1/canFam5 reference genome build). Additional genotyping confirmed a strong association between the 12-bp deletion and EOAD, where all affected dogs were homozygous for the deletion, while none of the control dogs was a deletion homozygote. A segregation pattern of this deletion in a 2-generation nuclear family indicated an autosomal recessive mode of inheritance. Since EPS8L2 plays a critical role in the maintenance and integrity of the inner ear hair cells in humans and other mammals, the inframe deletion found in this study represents a strong candidate causal mutation for EOAD in Rhodesian Ridgebacks. Genetic and clinical similarities between childhood deafness in humans and EOAD in Rhodesian Ridgebacks emphasizes the potential value of this dog breed in translational research in hereditary hearing disorders.

Genomic and Transcriptomic Characterization of Atypical Recurrent Flank Alopecia in the Cesky Fousek.

Non-inflammatory alopecia is a frequent skin problem in dogs, causing damaged coat integrity and compromised appearance of affected individuals. In this study, we examined the Cesky Fousek breed, which displays atypical recurrent flank alopecia (aRFA) at a high frequency. This type of alopecia can be quite severe and is characterized by seasonal episodes of well demarcated alopecic areas without hyperpigmentation. The genetic component responsible for aRFA remains unknown. Thus, here we aimed to identify variants involved in aRFA using a combination of histological, genomic, and transcriptomic data. We showed that aRFA is histologically similar to recurrent flank alopecia, characterized by a lack of anagen hair follicles and the presence of severely shortened telogen or kenogen hair follicles. We performed a genome-wide association study (GWAS) using 216 dogs phenotyped for aRFA and identified associations on chromosomes 19, 8, 30, 36, and 21, highlighting 144 candidate genes, which suggests a polygenic basis for aRFA. By comparing the skin cell transcription pattern of six aRFA and five control dogs, we identified 236 strongly differentially expressed genes (DEGs). We showed that the GWAS genes associated with aRFA are often predicted to interact with DEGs, suggesting their joint contribution to the development of the disease. Together, these genes affect four major metabolic pathways connected to aRFA: collagen formation, muscle structure/contraction, lipid metabolism, and the immune system.

Five genetic variants explain over 70% of hair coat pheomelanin intensity variation in purebred and mixed breed domestic dogs.

In mammals, the pigment molecule pheomelanin confers red and yellow color to hair, and the intensity of this coloration is caused by variation in the amount of pheomelanin. Domestic dogs exhibit a wide range of pheomelanin intensity, ranging from the white coat of the Samoyed to the deep red coat of the Irish Setter. While several genetic variants have been associated with specific coat intensity phenotypes in certain dog breeds, they do not explain the majority of phenotypic variation across breeds. In order to gain further insight into the extent of multigenicity and epistatic interactions underlying coat pheomelanin intensity in dogs, we leveraged a large dataset obtained via a direct-to-consumer canine genetic testing service. This consisted of genome-wide single nucleotide polymorphism (SNP) genotype data and owner-provided photos for 3,057 pheomelanic mixed breed and purebred dogs from 63 breeds and varieties spanning the full range of canine coat pheomelanin intensity. We first performed a genome-wide association study (GWAS) on 2,149 of these dogs to search for additional genetic variants that underlie intensity variation. GWAS identified five loci significantly associated with intensity, of which two (CFA15 29.8 Mb and CFA20 55.8 Mb) replicate previous findings and three (CFA2 74.7 Mb, CFA18 12.9 Mb, CFA21 10.9 Mb) have not previously been reported. In order to assess the combined predictive power of these loci across dog breeds, we used our GWAS data set to fit a linear model, which explained over 70% of variation in coat pheomelanin intensity in an independent validation dataset of 908 dogs. These results introduce three novel pheomelanin intensity loci, and further demonstrate the multigenic nature of coat pheomelanin intensity determination in domestic dogs.

Long-read assembly of a Great Dane genome highlights the contribution of GC-rich sequence and mobile elements to canine genomes.

Technological advances have allowed improvements in genome reference sequence assemblies. Here, we combined long- and short-read sequence resources to assemble the genome of a female Great Dane dog. This assembly has improved continuity compared to the existing Boxer-derived (CanFam3.1) reference genome. Annotation of the Great Dane assembly identified 22,182 protein-coding gene models and 7,049 long noncoding RNAs, including 49 protein-coding genes not present in the CanFam3.1 reference. The Great Dane assembly spans the majority of sequence gaps in the CanFam3.1 reference and illustrates that 2,151 gaps overlap the transcription start site of a predicted protein-coding gene. Moreover, a subset of the resolved gaps, which have an 80.95% median GC content, localize to transcription start sites and recombination hotspots more often than expected by chance, suggesting the stable canine recombinational landscape has shaped genome architecture. Alignment of the Great Dane and CanFam3.1 assemblies identified 16,834 deletions and 15,621 insertions, as well as 2,665 deletions and 3,493 insertions located on secondary contigs. These structural variants are dominated by retrotransposon insertion/deletion polymorphisms and include 16,221 dimorphic canine short interspersed elements (SINECs) and 1,121 dimorphic long interspersed element-1 sequences (LINE-1_Cfs). Analysis of sequences flanking the 3' end of LINE-1_Cfs (i.e., LINE-1_Cf 3'-transductions) suggests multiple retrotransposition-competent LINE-1_Cfs segregate among dog populations. Consistent with this conclusion, we demonstrate that a canine LINE-1_Cf element with intact open reading frames can retrotranspose its own RNA and that of a SINEC_Cf consensus sequence in cultured human cells, implicating ongoing retrotransposon activity as a driver of canine genetic variation.

R-locus for roaned coat is associated with a tandem duplication in an intronic region of USH2A in dogs and also contributes to Dalmatian spotting.

Structural variations (SVs) represent a large fraction of all genetic diversity, but how this genetic diversity is translated into phenotypic and organismal diversity is unclear. Explosive diversification of dog coat color and patterns after domestication can provide a unique opportunity to explore this question; however, the major obstacle is to efficiently collect a sufficient number of individuals with known phenotypes and genotypes of hundreds of thousands of markers. Using customer-provided information about coat color and patterns of dogs tested on a commercial canine genotyping platform, we identified a genomic region on chromosome 38 that is strongly associated with a mottled coat pattern (roaning) by genome-wide association study. We identified a putative causal variant in this region, an 11-kb tandem duplication (11,131,835-11,143,237) characterized by sequence read coverage and discordant reads of whole-genome sequence data, microarray probe intensity data, and a duplication-specific PCR assay. The tandem duplication is in an intronic region of usherin gene (USH2A), which was perfectly associated with roaning but absent in non-roaned dogs. We detected strong selection signals in this region characterized by reduced nucleotide diversity (π), increased runs of homozygosity, and extended haplotype homozygosity in Wirehaired Pointing Griffons and Australian Cattle Dogs (typically roaned breeds), as well as elevated genetic difference (FST) between Wirehaired Pointing Griffon (roaned) and Labrador Retriever (non-roaned). Surprisingly, all Dalmatians (N = 262) carried the duplication embedded in identical or similar haplotypes with roaned dogs, indicating this region as a shared target of selection during the breed's formation. We propose that the Dalmatian's unique spots were a derived coat pattern by establishing a novel epistatic interaction between roaning "R-locus" on chromosome 38 and an uncharacterized modifier locus. These results highlight the utility of consumer-oriented genotype and phenotype data in the discovery of genomic regions contributing to phenotypic diversity in dogs.

Examination of the efficacy of small genetic panels in genomic conservation of companion animal populations.

In many ways, dogs are an ideal model for the study of genetic erosion and population recovery, problems of major concern in the field of conservation genetics. Genetic diversity in many dog breeds has been declining systematically since the beginning of the 1800s, when modern breeding practices came into fashion. As such, inbreeding in domestic dog breeds is substantial and widespread and has led to an increase in recessive deleterious mutations of high effect as well as general inbreeding depression. Pedigrees can in theory be used to guide breeding decisions, though are often incomplete and do not reflect the full history of inbreeding. Small microsatellite panels are also used in some cases to choose mating pairs to produce litters with low levels of inbreeding. However, the long-term impact of such practices has not been thoroughly evaluated. Here, we use forward simulation on a model of the dog genome to examine the impact of using limited marker panels to guide pairwise mating decisions on genome-wide population-level genetic diversity. Our results suggest that in unmanaged populations, where breeding decisions are made at the pairwise-rather than population-level, such panels can lead to accelerated loss of genetic diversity at genome regions unlinked to panel markers, compared to random mating. These results demonstrate the importance of genome-wide genetic panels for managing and conserving genetic diversity in dogs and other companion animals.

Body size, inbreeding, and lifespan in domestic dogs.

Inbreeding poses a real or potential threat to nearly every species of conservation concern. Inbreeding leads to loss of diversity at the individual level, which can cause inbreeding depression, and at the population level, which can hinder ability to respond to a changing environment. In closed populations such as endangered species and ex situ breeding programs, some degree of inbreeding is inevitable. It is therefore vital to understand how different patterns of breeding and inbreeding can affect fitness in real animals. Domestic dogs provide an excellent model, showing dramatic variation in degree of inbreeding and in lifespan, an important aspect of fitness that is known to be impacted by inbreeding in other species. There is a strong negative correlation between body size and lifespan in dogs, but it is unknown whether the higher rate of aging in large dogs is due to body size per se or some other factor associated with large size. We used dense genome-wide SNP array data to calculate average inbreeding for over 100 dog breeds based on autozygous segment length and found that large breeds tend to have higher coefficients of inbreeding than small breeds. We then used data from the Veterinary medical Database and other published sources to estimate life expectancies for pure and mixed breed dogs. When controlling for size, variation in inbreeding was not associated with life expectancy across breeds. When comparing mixed versus purebred dogs, however, mixed breed dogs lived about 1.2 years longer on average than size-matched purebred dogs. Furthermore, individual pedigree coefficients of inbreeding and lifespans for over 9000 golden retrievers showed that inbreeding does negatively impact lifespan at the individual level. Registration data from the American Kennel Club suggest that the molecular inbreeding patterns observed in purebred dogs result from specific breeding practices and/or founder effects and not the current population size. Our results suggest that recent inbreeding, as reflected in variation within a breed, is more likely to affect fitness than historic inbreeding, as reflected in variation among breeds. Our results also indicate that occasional outcrosses, as in mixed breed dogs, can have a substantial positive effect on fitness.

A genome-wide association study of deafness in three canine breeds.

Congenital deafness in the domestic dog is usually related to the presence of white pigmentation, which is controlled primarily by the piebald locus on chromosome 20 and also by merle on chromosome 10. Pigment-associated deafness is also seen in other species, including cats, mice, sheep, alpacas, horses, cows, pigs, and humans, but the genetic factors determining why some piebald or merle dogs develop deafness while others do not have yet to be determined. Here we perform a genome-wide association study (GWAS) to identify regions of the canine genome significantly associated with deafness in three dog breeds carrying piebald: Dalmatian, Australian cattle dog, and English setter. We include bilaterally deaf, unilaterally deaf, and matched control dogs from the same litter, phenotyped using the brainstem auditory evoked response (BAER) hearing test. Principal component analysis showed that we have different distributions of cases and controls in genetically distinct Dalmatian populations, therefore GWAS was performed separately for North American and UK samples. We identified one genome-wide significant association and 14 suggestive (chromosome-wide) associations using the GWAS design of bilaterally deaf vs. control Australian cattle dogs. However, these associations were not located on the same chromosome as the piebald locus, indicating the complexity of the genetics underlying this disease in the domestic dog. Because of this apparent complex genetic architecture, larger sample sizes may be needed to detect the genetic loci modulating risk in piebald dogs.

Genetic mapping of distal femoral, stifle, and tibial radiographic morphology in dogs with cranial cruciate ligament disease.

Cranial cruciate ligament disease (CCLD) is a complex trait. Ten measurements were made on orthogonal distal pelvic limb radiographs of 161 pure and mixed breed dogs with, and 55 without, cranial cruciate partial or complete ligament rupture. Dogs with CCLD had significantly smaller infrapatellar fat pad width, higher average tibial plateau angle, and were heavier than control dogs. The first PC weightings captured the overall size of the dog's stifle and PC2 weightings reflected an increasing tibial plateau angle coupled with a smaller fat pad width. Of these dogs, 175 were genotyped, and 144,509 polymorphisms were used in a genome-wide association study with both a mixed linear and a multi-locus model. For both models, significant (pgenome <3.46×10-7 for the mixed and< 6.9x10-8 for the multilocus model) associations were found for PC1, tibial diaphyseal length and width, fat pad base length, and femoral and tibial condyle width at LCORL, a known body size-regulating locus. Other body size loci with significant associations were growth hormone 1 (GH1), which was associated with the length of the fat pad base and the width of the tibial diaphysis, and a region on CFAX near IRS4 and ACSL4 in the multilocus model. The tibial plateau angle was associated significantly with a locus on CFA10 in the linear mixed model with nearest candidate genes BET1 and MYH9 and on CFA08 near candidate genes WDHD1 and GCH1. MYH9 has a major role in osteoclastogenesis. Our study indicated that tibial plateau slope is associated with CCLD and a compressed infrapatellar fat pad, a surrogate for stifle osteoarthritis. Because of the association between tibial plateau slope and CCLD, and pending independent validation, these candidate genes for tibial plateau slope may be tested in breeds susceptible to CCLD before they develop disease or are bred.

Imputation of canine genotype array data using 365 whole-genome sequences improves power of genome-wide association studies.

Genomic resources for the domestic dog have improved with the widespread adoption of a 173k SNP array platform and updated reference genome. SNP arrays of this density are sufficient for detecting genetic associations within breeds but are underpowered for finding associations across multiple breeds or in mixed-breed dogs, where linkage disequilibrium rapidly decays between markers, even though such studies would hold particular promise for mapping complex diseases and traits. Here we introduce an imputation reference panel, consisting of 365 diverse, whole-genome sequenced dogs and wolves, which increases the number of markers that can be queried in genome-wide association studies approximately 130-fold. Using previously genotyped dogs, we show the utility of this reference panel in identifying potentially novel associations, including a locus on CFA20 significantly associated with cranial cruciate ligament disease, and fine-mapping for canine body size and blood phenotypes, even when causal loci are not in strong linkage disequilibrium with any single array marker. This reference panel resource will improve future genome-wide association studies for canine complex diseases and other phenotypes.

Inbreeding depression causes reduced fecundity in Golden Retrievers.

Inbreeding depression has been demonstrated to impact vital rates, productivity, and performance in human populations, wild and endangered species, and in recent years, the domestic species. In all cases, standardized, high-quality phenotype data on all individuals are invaluable for longitudinal analyses such as those required to evaluate vital rates of a study cohort. Further, many investigators agree upon the preference for and utility of genomic measures of inbreeding in lieu of pedigree-based estimates of inbreeding. We evaluated the association of measures of reproductive fitness in 93 Golden Retrievers enrolled in the Golden Retriever Lifetime Study with a genomic measurement of inbreeding, FROH. We demonstrate a statistically significant negative correlation between fecundity and FROH. This work sets the stage for larger scale analyses to investigate genomic regions associated with fecundity and other measures of fitness.

Cellular energetics and mitochondrial uncoupling in canine aging.

The first domesticated companion animal, the dog, is currently represented by over 190 unique breeds. Across these numerous breeds, dogs have exceptional variation in lifespan (inversely correlated with body size), presenting an opportunity to discover longevity-determining traits. We performed a genome-wide association study on 4169 canines representing 110 breeds and identified novel candidate regulators of longevity. Interestingly, known functions within the identified genes included control of coat phenotypes such as hair length, as well as mitochondrial properties, suggesting that thermoregulation and mitochondrial bioenergetics play a role in lifespan variation. Using primary dermal fibroblasts, we investigated mitochondrial properties of short-lived (large) and long-lived (small) dog breeds. We found that cells from long-lived breeds have more uncoupled mitochondria, less electron escape, greater respiration, and capacity for respiration. Moreover, our data suggest that long-lived breeds have higher rates of catabolism and ß-oxidation, likely to meet elevated respiration and electron demand of their uncoupled mitochondria. Conversely, cells of short-lived (large) breeds may accumulate amino acids and fatty acid derivatives, which are likely used for biosynthesis and growth. We hypothesize that the uncoupled metabolic profile of long-lived breeds likely stems from their smaller size, reduced volume-to-surface area ratio, and therefore a greater need for thermogenesis. The uncoupled energetics of long-lived breeds lowers reactive oxygen species levels, promotes cellular stress tolerance, and may even prevent stiffening of the actin cytoskeleton. We propose that these cellular characteristics delay tissue dysfunction, disease, and death in long-lived dog breeds, contributing to canine aging diversity.

Fine-Scale Resolution of Runs of Homozygosity Reveal Patterns of Inbreeding and Substantial Overlap with Recessive Disease Genotypes in Domestic Dogs.

Inbreeding leaves distinct genomic traces, most notably long genomic tracts that are identical by descent and completely homozygous. These runs of homozygosity (ROH) can contribute to inbreeding depression if they contain deleterious variants that are fully or partially recessive. Several lines of evidence have been used to show that long (> 5 megabase) ROH are disproportionately likely to harbor deleterious variation, but the extent to which long vs. short tracts contribute to autozygosity at loci known to be deleterious and recessive has not been studied. In domestic dogs, nearly 200 mutations are known to cause recessive diseases, most of which can be efficiently assayed using SNP arrays. By examining genome-wide data from over 200,000 markers, including 150 recessive disease variants, we built high-resolution ROH density maps for nearly 2,500 dogs, recording ROH down to 500 kilobases. We observed over 678 homozygous deleterious recessive genotypes in the panel across 29 loci, 90% of which overlapped with ROH inferred by GERMLINE. Although most of these genotypes were contained in ROH over 5 Mb in length, 14% were contained in short (0.5 - 2.5 megabase) tracts, a significant enrichment compared to the genetic background, suggesting that even short tracts are useful for computing inbreeding metrics like the coefficient of inbreeding estimated from ROH (FROH ). In our dataset, FROH differed significantly both within and among dog breeds. All breeds harbored some regions of reduced genetic diversity due to drift or selective sweeps, but the degree of inbreeding and the proportion of inbreeding caused by short vs. long tracts differed between breeds, reflecting their different population histories. Although only available for a few species, large genome-wide datasets including recessive disease variants hold particular promise not only for disentangling the genetic architecture of inbreeding depression, but also evaluating and improving upon current approaches for detecting ROH.

Direct-to-consumer DNA testing of 6,000 dogs reveals 98.6-kb duplication associated with blue eyes and heterochromia in Siberian Huskies.

Consumer genomics enables genetic discovery on an unprecedented scale by linking very large databases of personal genomic data with phenotype information voluntarily submitted via web-based surveys. These databases are having a transformative effect on human genomics research, yielding insights on increasingly complex traits, behaviors, and disease by including many thousands of individuals in genome-wide association studies (GWAS). The promise of consumer genomic data is not limited to human research, however. Genomic tools for dogs are readily available, with hundreds of causal Mendelian variants already characterized, because selection and breeding have led to dramatic phenotypic diversity underlain by a simple genetic structure. Here, we report the results of the first consumer genomics study ever conducted in a non-human model: a GWAS of blue eyes based on more than 3,000 customer dogs with validation panels including nearly 3,000 more, the largest canine GWAS to date. We discovered a novel association with blue eyes on chromosome 18 (P = 1.3x10-68) and used both sequence coverage and microarray probe intensity data to identify the putative causal variant: a 98.6-kb duplication directly upstream of the Homeobox gene ALX4, which plays an important role in mammalian eye development. This duplication is largely restricted to Siberian Huskies, is strongly associated with the blue-eyed phenotype (chi-square P = 5.2x10-290), and is highly, but not completely, penetrant. These results underscore the power of consumer-data-driven discovery in non-human species, especially dogs, where there is intense owner interest in the personal genomic information of their pets, a high level of engagement with web-based surveys, and an underlying genetic architecture ideal for mapping studies.

The evolutionary history of dogs in the Americas.

Dogs were present in the Americas before the arrival of European colonists, but the origin and fate of these precontact dogs are largely unknown. We sequenced 71 mitochondrial and 7 nuclear genomes from ancient North American and Siberian dogs from time frames spanning ~9000 years. Our analysis indicates that American dogs were not derived from North American wolves. Instead, American dogs form a monophyletic lineage that likely originated in Siberia and dispersed into the Americas alongside people. After the arrival of Europeans, native American dogs almost completely disappeared, leaving a minimal genetic legacy in modern dog populations. The closest detectable extant lineage to precontact American dogs is the canine transmissible venereal tumor, a contagious cancer clone derived from an individual dog that lived up to 8000 years ago.

Elucidating biogeographical patterns in Australian native canids using genome wide SNPs.

Dingoes play a strong role in Australia's ecological framework as the apex predator but are under threat from hybridization and agricultural control programs. Government legislation lists the conservation of the dingo as an important aim, yet little is known about the biogeography of this enigmatic canine, making conservation difficult. Mitochondrial and Y chromosome DNA studies show evidence of population structure within the dingo. Here, we present the data from Illumina HD canine chip genotyping for 23 dingoes from five regional populations, and five New Guinea Singing Dogs to further explore patterns of biogeography using genome-wide data. Whole genome single nucleotide polymorphism (SNP) data supported the presence of three distinct dingo populations (or ESUs) subject to geographical subdivision: southeastern (SE), Fraser Island (FI) and northwestern (NW). These ESUs should be managed discretely. The FI dingoes are a known reservoir of pure, genetically distinct dingoes. Elevated inbreeding coefficients identified here suggest this population may be genetically compromised and in need of rescue; current lethal management strategies that do not consider genetic information should be suspended until further data can be gathered. D statistics identify evidence of historical admixture or ancestry sharing between southeastern dingoes and South East Asian village dogs. Conservation efforts on mainland Australia should focus on the SE dingo population that is under pressure from domestic dog hybridization and high levels of lethal control. Further data concerning the genetic health, demographics and prevalence of hybridization in the SE and FI dingo populations is urgently needed to develop evidence based conservation and management strategies.

Comparison of village dog and wolf genomes highlights the role of the neural crest in dog domestication.

BACKGROUND: Domesticated from gray wolves between 10 and 40 kya in Eurasia, dogs display a vast array of phenotypes that differ from their ancestors, yet mirror other domesticated animal species, a phenomenon known as the domestication syndrome. Here, we use signatures persisting in dog genomes to identify genes and pathways possibly altered by the selective pressures of domestication. RESULTS: Whole-genome SNP analyses of 43 globally distributed village dogs and 10 wolves differentiated signatures resulting from domestication rather than breed formation. We identified 246 candidate domestication regions containing 10.8 Mb of genome sequence and 429 genes. The regions share haplotypes with ancient dogs, suggesting that the detected signals are not the result of recent selection. Gene enrichments highlight numerous genes linked to neural crest and central nervous system development as well as neurological function. Read depth analysis suggests that copy number variation played a minor role in dog domestication. CONCLUSIONS: Our results identify genes that act early in embryogenesis and can confer phenotypes distinguishing domesticated dogs from wolves, such as tameness, smaller jaws, floppy ears, and diminished craniofacial development as the targets of selection during domestication. These differences reflect the phenotypes of the domestication syndrome, which can be explained by alterations in the migration or activity of neural crest cells during development. We propose that initial selection during early dog domestication was for behavior, a trait influenced by genes which act in the neural crest, which secondarily gave rise to the phenotypes of modern dogs.